Exact Electron-Pairing Ground States of Tight-Binding Models with Local Attractive 

Interactions 
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We present a class of exactly solvable models of correlated electrons. The models are defined 
in any dimension d and consist of electron-hopping terms and local attractive interactions between 
electrons. For each even number of electrons less than or equal to l/(d + l)-filling, we find the exact 
ground state in which all electrons form pairs of a certain type, and thus the models exhibit an 
electron-pair condensation. 

PACS numbers: 71.10.Fd, 74.20.-z 



The origin of high-temperature superconductivity has 
attracted much interest and is still controversial. On the 
other hand, the properties of usual (low-temperature) 
superconductors are well explained by the celebrated 
Bardeen-Cooper-Schrieffcr (BCS) theory l]. On the ba- 
sis of the BCS theory, one treats electron systems with ef- 
fective electron-electron attraction arising from electron- 
phonon couplings and discusses a phase transition asso- 
ciated with an electron-pair condensation by using mean- 
field type approximations. However, whether attractive 
interactions between electrons really induce the forma- 
tion of many electron pairs and their condensation is not 
at all trivial, and it is desirable to clarify this point with- 
out relying on approximation methods. 

One of the results in this direction was obtained by 
Shen and Qiu Q , who proved the existence of off-diagonal 
long-range order, closely related to a condensation, in 
the Hubbard model with on-site attractive interactions 
on certain bipartite lattices. Another was obtained by 
Essler, Korepin, and Schoutens 0], who proposed solv- 
able extended Hubbard models exhibiting a condensation 
of electrons forming the r\ pair j^, |5j which is a superpo- 
sition of on-site pairs with momentum 0. There are some 
rigorous results related to these models, in particular, in 
one dimension where the models are exactly solvable by 
Bethe ansatz (see, for example, 0, S an d references 
therein). Results in higher dimensions are limited. 

In this letter we introduce a new class of exactly solv- 
able models exhibiting a condensation at zero tempera- 
ture 0, llfll | . The models are defined in arbitrary dimen- 
sions and consist of hopping terms of electrons and local 
attractive interactions between electrons. We show that 
the models have the ground states in which all electrons 
form pairs of a certain type. It is worth to note that 
our pair, unlike the r\ pair used in the results mentioned 
above, is made up of spin-singlet pairs of electrons at 
different sites in addition to on-site pairs. We also show 
that, in the momentum space, the exact ground states are 
represented as superpositions of products of pairs formed 
from electrons with opposite wave vectors and spins. 

As seen later, the present models are related to nearly- 
flat-band models, which are proved to exhibit ferromag- 



netism Q 0. Several authors d H E3 
proposed possibilities of experimental realization of the 
nearly-flat-band ferromagnetism. Our results suggest 
that there is a chance of finding superconductivity or 
interesting phenomena caused by competition between 
magnetism and superconductivity as well as ferromag- 
netism in the proposed systems. 

To simplify the discussion, we describe the results in 
the simplest version of the models. We can consider 
similar models on other lattices constructed by using 
the same methods as that by Mielke ^lj and that by 
Tasaki , which were developed in the study of (nearly) 
flat-band ferromagnetism. Let V be a (i-dimcnsional hy- 
percubic lattice [1, L] d n Z d with periodic boundary con- 
ditions, where L is an arbitrary integer. Let M be a 
collection of sites located at the mid-points of nearest- 
neighbor pairs in V. Then we define decorated lattice 
A = V U M. We denote by Cj iCT and cj a the annihilation 
and the creation operators, respectively, for an electron 
with spin <j =\,\ at site i G A. These operators satisfy 
usual fermion anticommutation relations. The number 
operator m, a is defined as — c} ia Ci. a . We denote 
by N c the electron number and by <&o a state without 
electrons. We define fermion operators, which play an 
essential role in our model, as 



«£iW;|«-x|=l/2 



for x £ V and 
6 



II. (7 ~ C U,(T H~ OL ^ C'X,<7 

xeV;\x-u\=l/2 



(i) 



(2) 



for u S M, where a is a real number. We note that lo- 
calized single-electron states corresponding to a-operator 
and 6-operator are orthogonal because 



{b Ut(r ,al } = 



(3) 



for any x S V and any u G M . 

We consider the tight-binding model of electrons with 
local attractive interactions described by the Hamilto- 
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(a) 




(b) 




FIG. 1: The model in one dimension, (a) The lattice structure 
and hopping matrix elements. The solid and open circles 
represent sites in V and M, respectively. The on-site potential 
tu is s + 1c? t if i € V and 2a 2 s + t if z e M. The hopping 

a 2 t if \i - .7 1 = l,i, 7 G V, 
J\ = 



amplitudes are given by ti 



tij = a s if \i — j\ = 1, i, j £ M, tij = —a(s 
1/2, and tij = otherwise 



= 1, *% J 
t) if |t - 



(b) The dispersion relations. 
The operator aJ. j(T creates an electron in the band with the 
dispersion relation £ a (k) = s + 2q 2 s(1 + cos k) and b„ :<T creates 
an electron in the band with et,i(k) = t + 2a 2 t(l + cos k). 



nian H = H£ op 



H\„t with 



^hop — S a x,o- a ' x >< T > 
xEV<r=1,l 

< P = t E E b l* b »,<" 

Hint = -W ^ al^al d a Xtl a Xt 

x£V 



(4) 
(5) 
(6) 



where s,t and W are positive parameters. It is noted 
that H conserves the electron number N e and possesses 
the spin SU(2) symmetry. 

The sum of H£ op and H^ op can be reduced to a stan- 



dard tight-binding Hamiltonian ^2% ,gA So- * 



describing quantum mechanical motion of electrons in 
A. It follows from that the dispersion relations of 
the model are classified into two types; one corresponds 
to single-electron states {a x ^o)x£V , and the others cor- 
respond to (r $o)'ueM- More precisely, the dispersion 
relations are given by e a {k) = s + 2a 2 s 5^ J _ 1 (1 + cos ki), 
£bs(k) = t + 2a 2 t X)f=i(l + coski), and £& >m (fe) — t with 
m = 2, . . . , d (see Fig. 1 for an example in one dimen- 
sion). Here, k is the wave vector in the set 0] 

K = {(27TTU/L, . . . , 2Tm d /L)\m = 0, ±1, . . . , ±{L - l)/2} . 

(7) 

The energy eigenstate with eigenvalue s a {k) is given by 



L*o with 



£a(k)\V\ 



E eife 'H, 



(8) 



The many-electron ground state of is usually 

the Fermi sea [2p| . 

By rewriting each term in iJj nt as 

-W \\a\\ 4 + W \\af a x , T 4 )t + H^a^a^a^ (9) 

with ||o|| = (l + 2da 2 ), one finds that this is bounded be- 
low by — W \ \a\\ , which is attained by states in the form 
of a x ^a x I • • • $o • The terms in H mi are thus interpreted 
as attractive interactions between electrons in localized 
single-electron states corresponding to a-operator. In the 
/c-space, Hint is represented as 



1 

rn 



E W k,k'A a l+ q ,\ a k'-q,i a k',i a Kh ( 10 ) 



k,k',qelC 
W 



W kjk ^ q = -We a (/c + q)e a {k' - q)e a {k')e a (k), (11) 

which has a similar form to the BCS interaction Hamilto- 
nian |2jj , although there is a difference in the dependence 
of the amplitudes on the wave vectors. It is noted that 
interaction H mt contains scattering processes of electron 
pairs with non-zero total momentum, which are usually 
dropped in the BCS theory. Nevertheless, our theorem 
below establishes that the ground state of H with fixed 
even N Q consists of electron pairs having total momen- 
tum (see JTHJ and fity ). 

Before stating our main result, we need to introduce 
further notation. Let G be the \V\ x \V\ Gram ma- 
trix for a-operator, whose matrix elements are defined 
by (G)xy = {a x , CT , a\ a } for x,y e V. Since (a x<T $o)xeV 
is linearly independent, G is regular and has its inverse 
matrix G _1 . Then we introduce the dual operators of o- 
operator as a Xj(T = Ylyevi^ ) x v a y> rT - ^ * s eas Y to check 
the inverse relation a x a — X^yey {G) xy a y ^ a and the anti- 
commutation relations 



for x, y € V and a, t =t, j. Let us define 



/ 4 (^)^y^x,t%,i^ 

x,y£V 



(12) 



(13) 



which creates two electrons in a spin-singlet pairing state. 
We note that is rewritten as 



c 



1 - E h l,A,i 



E a i,T a i,i 

x£V 



16V x£V x,y£V 

(14) 

since G and G 1 are real symmetric matrices. Our main 
result in this letter is summarized as follows: 
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Theorem. Consider Hamiltonian H with W — 2s/||a|| 
and fixed electron number N c in N c < 2\V\. For even 
N c , the ground state &g.n is unique, has zero energy, 
and is expressed as 



$G,iV e = (C f ) * *0- 



(15) 



When N c is odd, the ground state has positive energy. 

It is noted that the ground state energy for odd N e 
may converge to zero as L — > oo. Whether this is the 
case or not should be clarified in a future study. 

If we set s < and W = 0, Hamiltonian H = 
H£ op + -f/^op i s ec l ua l to the hop ping term of the nearly- 
flat-band model studied in Ref. [12J. Thus, in this case, 
the model has saturated ferromagnetic ground states for 
N e = \V\ when sufficiently large on-site repulsion is 
added. We note that the ferromagnetic ground states 
found in the above situation survive for all W < and 
also for sufficiently small positive values of W [22j • The 
present model with the on-site repulsion is expected to 
describe a quantum phase transition between magnetic 
states and electron-pairing states. 

The elements of Gram matrix G are explicitly given by 
(G) xy = 1 + 2da 2 ifx = y, (G) xy = a 2 if \x - y\ = 1, 
and (G) xy = otherwise, and a standard Fourier analysis 
yields (G-% y = \V\~ l j: keK7 ^e ik <*-y\ From this, 

we find that the pair creation operator is written as 



(16) 



Substituting this expression into (|15|) . one finds that 
the ground state is a superposition of products of pairs 
a \ T a -fe J." T° see the representation of C} in terms of 
c-operator, we substitute a\. a — X^eA^^ ! a m t° the 
final expression in 114J1 . where ip X i = {a Xi(J , cj a }. Then 
we obtain 



c 



t — \ " t t 

ieA 



E 



(17) 



with Wij = Y,xjjev(G 1 )xy<Pxi<P y j- Each term in the 
second sum in (|17|) corresponds to a spin-singlet pair of 
electrons at different sites. In the case of x, y € V, the 
weight w xy of a pair reduces to (G -1 )^, which has the 
Ornstein-Zernicke behavior for \x — y\ 1. 

By making a linear combination of &G.N a with different 
electron numbers, we can construct an explicit electron- 
number symmetry breaking ground state. To see this, let 
us introduce order parameter Aa = £/|V|. A straight- 
forward calculation yields 



A.N C 



($ GiArc ,A A A A $ G 

($G,JV ,$G,7V,) 



fJ-A,N c (18) 



with/x AiiVe = (N c /(2\V\) + l/\V\)(l- N C /(2\V)). Taking 
a limit L, N c — > oo so that the electron density N c / (2|A|) 



will converge to v, we conclude that ( AaAI ) 



with fj,„ = {d+l)v{l— (d+l)v). Then, we define $ G Wo = 

$G jv , which is also a zero- 



A, At) 



1 + a a/\J\ la a^ A/ 
energy state of H . It is easy to see that 

(*G,JV.. A A*G,W.) 



(A 



A /A,N B 



( $ G,JV 6 > $ G,JvJ 



W2, (19) 



which implies the explicit electron-number symmetry 
breaking for < v < l/(d + 1). 

Proof of the theorem. By using W — 2s/||a|| 2 , we 
rewrite H as H = H^ op + H' int with 



"int — 2 Z j / t x - a 

X £V <T=1d 



X, — <T Uj X ^ — (T Uj X.(T- 



(20) 



(Here — <r denotes spin opposite to a.) This expression 
is crucial in our proof. Since both H^ op an< ^ ^int are 
positive semidefinite operators, one can conclude that a 
zero-energy state of both of these terms is a ground state. 
We first show that $g,jv c m lfT3|) is in fact a zero-energy 
state of H. From lfP2)l and {a\. a ) 2 = 0, one finds that 

4,4. a *,tC t = C + 4,4. a *,t + ai,i a i,i = a l,l a ^,7- ( 21 ) 

Similarly, we have a} x *a x ,iC^ — C^ a l ^ a x.i- These two re- 
lations imply that iJ int $G,Af = 0- Furthermore it imme- 
diately follows from © that ff£ $ GiJ v c = 0. Therefore 
<&G,jv is a ground state having zero energy. 

To prove the other statement in the theorem, we use 
the following lemma, which will be proved later. 
Lemma. Any zero-energy state <I> of H with N c < 2\V\ 
(where N c is not fixed) is written as 




(22) 



AcV 



where coefficients 4>a satisfy 4>a — <pA' if \A\ = l-^'l- 
This lemma implies that the ground state energy for fixed 
odd N c is positive. Let us prove the uniqueness of the 
ground state for fixed even N c . Now suppose that there 
are two zero-energy states. Then any linear combination 
of these states is also a zero-energy state, which must 
satisfy the statement in the lemma. However, this implies 
that all the coefficients are vanishing, since we can make 
a suitable linear combination such that a coefficient <f>A 
for a subset Aq is zero. This is contradicting with the 
assumption, and thus the ground state, which must have 
zero energy, is unique. This completes the proof of the 
theorem. ■ 

Proof of the lemma. Let $ be an arbitrary zero-energy 
state of H with N c < 2\V\. Since each term tb' u a b u . a in 
-ffhop ^ s positive semidefinite, $ should satisfy = 



4 



for any u £ M and a =f , J.. Thus $ must be of the form does not contain a site y with |y 

J. 

(A\{x},AU{</})' 



1. By checking 



the coefficient of $ 



we have 



e n n 



$0, 



A T ,AiCV 



(23) 

where $(^-.4,) are suitable coefficients. 

To be a zero-energy state, $ must further- 
more satisfy H' int & = 0, i.e., a x _ a a x ^<& = 
for any x £ V and a = f , |. By using 

*(A t ,A 4 ) = (fLeA, 4,t) (iLeA^ 4,J $ o, we expand $ 



At.A^cV 



(24) 



^(sgnfy, x; A]cf> A + sgn[x, y; A x . 



-yi 



.,) = 0, (29) 

where A K ^ y = (A\{x}) U {?/}. Since F* x for 
|x — y\ = 1 is non-zero by definition and sgn[y, x; A] — 
— sgnfa;, y; A x ^ y ], we find that the coefficients must sat- 
isfy 4>a = 4>A X ^ ■ Repeating the same argument for all 
x € V, we conclude that <j> A — <f>A> whenever \A\ = \A'\. 
This completes the proof of the lemma. ■ 

I would like to thank Masanori Yamanaka for useful 
discussions. 



where 4>^ Ar A L ) are new coefficients, and operate a x ^a Xj 
on $ in this form. Then, by using (|12|l . we have 

X[x 6 A h x ^A i ]8ga[x;A r ,A l ] 



A T ,A±CV 



X(j) 



(A r ,Al)"(A r \{x},AlU{x}) 



= 0, 



(25) 



where sgn[- • • ] is a sign factor arising from exchanges of 
fermion operators, and x[" even t"] takes 1 if "event" is 
true and takes otherwise. Since all the terms in the left 
hand side in the above equation are linearly independent, 
we obtain (j>^ A A ) = if both x £ Af and x ^ Ap Since 
this must hold for all x £ V, we find that the sum in i|24|) 
is restricted to subsets Aj, Aj such that A-f C Ap 

By using ®{ AuAl) = (EUa t 4,t) (iLeA; 4,J $0 ' 
and by taking account of the above result, we again 
rewrite as 



2J 0(A T) Ai)*(A T ,Ai) 
A r ,A l CV;\A- t \<\A l \ 



(26) 



with new coefficients 



11 



"(a t ,a,)- °P eratin g a x,l a *,l on * 
in expression i|2ti|) and repeating a similar argument to 

the above, we find that cf>'^ A ^ A % must be zero if A-\ ^ Aj_. 
Thus we have shown that $ must be written as 



Z— < -- V (A,A)' 



ACV 



(27) 



where 0a = 4>( AA y Now we again operate a x ^a x ,i — 

Sy.y'ey(C)a:y( C? )^'4,i^''T 011 * in expression and 
derive conditions on <^4. The resulting equation is 



E E x[y^A,y'£A}sgn[y,y';A] 

ACV y,y'GV 



I 

(A\{y'},AU{y}) 



= 0, 



(28) 



where Fy y , = (G) xy (G) xy i and sgn[- • ■] is a fermion-sign 
factor. Let us choose a subset A which contains site x and 
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We assume that L is an odd integer here. 
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Note that (aj. ^o)keK is orthonormal. 
When s < 0, the model with N c — | V\ is proved to exhibit 
saturated ferromagnetism for sufficiently large negative 
values of W even in the absence of on-site repulsion. 



